DTP/00/34 
May 2000 



The unintegrated gluon distribution 
from the CCFM equation 



M.A. Kimber a , J. Kwiecinski a > fe , A.D. Martin a and A.M. Stasto a ' b 

a Department of Physics, University of Durham, Durham, DH1 3LE 
b H. Niewodniczanski Institute of Nuclear Physics, ul. Radzikowskiego 152, Krakow, Poland 



Abstract 

The gluon distribution f(x,k 2 ,/i 2 ), unintegrated over the transverse momentum kt of 
the gluon, satisfies the angular-ordered CCFM equation which interlocks the dependence 
on the scale kt with the scale \i of the probe. We show how, to leading log accuracy, the 
equation can be simplified to a single scale problem. In particular we demonstrate how 
to determine the two-scale unintegrated distribution f(x,k 2 ,fi 2 ) from knowledge of the 
integrated gluon obtained from a unified scheme embodying both BFKL (log(l/a;)) and 
DGLAP (log// 2 ) evolution. 



1 Introduction 



Deep inelastic electron-proton scattering is described in terms of scale dependent parton dis- 
tributions q(x,fi 2 ) and g(x,fi 2 ). For less inclusive processes it is however necessary to consider 
distributions unintegrated over the transverse momentum k t of the parton, which for the gluon, 
for example, we denote by f(x, k 2 ,fi 2 ). These distributions depend on two hard scales; k t and 
the hard scale \x of the probe. The (conventional) integrated gluon distribution is given by 

rfi 2 ^2 

xg(x,fi 2 ) = J f(x,k 2 ,fi 2 ). (1) 

Unintegrated distributions are required to describe measurements where transverse momenta 
are exposed explicitly. For example to describe the pt spectrum of prompt photons produced 
in high energy hadron collisions or for dijets or vector mesons produced at HERA. 
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At very low x, that is to leading log(l/x) accuracy, the unintegrated distribution becomes 
independent of the hard scale and so from ([!]) we have 



f(x, k 2 , [i 2 



d 



d In A 2 



(xg(x, A 2 ) 



(2) 



Clearly cannot remain true as x increases. Indeed we see that it would give negative values 
for /. Moreover even at low x, there are significant subleading corrections which, to some level 
of approximation, modify (@) to the form |], § 
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A 2 ) T g (\,n) 



(3) 



where T 9 is the Sudakov form factor. In fact (|3|) is oversimplified. As discussed below, the 
expression for f(x, k 2 ,fi 2 ) is more complicated than fl3|). 

The natural framework for unifying the small and large x domains is the CCFM formalism 
based on angular ordering Pfl-|JL0[|, which follows from colour coherence effects [IT]. It reduces 
to the leading order DGLAP formalism at moderate x and it embodies the BFKL formalism 
at small x. The unintegrated gluon distribution f(x, k 2 ,fi 2 ) satisfies the CCFM equation [|, [|] 
which interlocks the two hard scales (k 2 , /x 2 ) in a complicated way. The equation is based on the 
coherent radiation of gluons, which leads to an angular ordering of the gluon emissions along 
the chain. The ordering introduces a scale specifying the maximum angle of gluon emission, 
which turns out to be essentially the hard scale \i of the probe. At moderate x the angular 
ordering becomes an ordering in the gluon transverse momenta and the CCFM equation reduces 
to DGLAP evolution. At very small x the angular ordering does not provide any constraint on 
the transverse momenta along the chain and, in the leading log(l/a;) approximation, f(x, k 2 , /j, 2 ) 
becomes the /i independent distribution which satisfies the BFKL equation. On the other hand, 
although the dependence on the scale \x only enters at subleading log(l/x) level, / does depend 
on \i through leading log /i 2 evolution. 

The outline of the paper is as follows. The angular-ordered CCFM equation is introduced 
in Section 2. In Section 3 we simplify this evolution, yet staying within leading log accuracy, 
to show that the two-scale distribution f(x,k 2 ,fj, 2 ) can be obtained in terms of the conven- 
tional one-scale g(x,fi 2 ) distribution. In this way, we are led to a procedure for determining 
f(x, k 2 , ii 2 ) from a unified DGLAP/BFKL type single-scale evolution equation. This is de- 
scribed in Section 4. Moreover we are able to extend the formalism to incorporate important 
subleading log(l/x) effects, which are generated by the so-called consistency conditionF] |12|, [L3[ 
and which subsume the angular ordering constraint at low x. We also extend the formalism 
to include the contributions due to the quark distributions. For comparison, in Section 5 we 
present the pure DGLAP-type approach to determine f(x, k 2 ,fi 2 ), in which the gluon cascade 
evolves according to evolution strongly ordered in kt- Section 6 contains sample numerical 



Called the kinematic constraint in [fL2l 
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results for f(x, k 2 ,fi 2 ) obtained from the fully unified approach of Section 4. As expected, we 
have some diffusion of the gluon transverse momenta into the region k t > //. This is in contrast 
to the pure DGLAP-type approximation in which the distribution f(x, k 2 ,fi 2 ) is limited to the 
domain k t < fi. Finally, Section 7 contains a summary of the procedure that we have used to 
determine the unintegrated gluon distribution f(x, k 2 , /i 2 ). 

2 The CCFM equation 

The unintegrated gluon distribution satisfies the CCFM evolution equation || f|] based on 
angular ordering. In unfoldedQ form the equation is 

f(x,k 2 ,fi 2 ) = f (x,k 2 ) 

Q(z — x) 0(/i — qz) 

1 2 

|tp(z)/(£ *f,, 2 ) - zP(z)e(p-q)l(x,klq>) 

- sM't/S^-^-M!-* 9 ')- < 4 > 

where k' t = \k t + (1 — z)q\. For simplicity, as is taken outside the integrals, but the scale will 
be specified carefully in the final equations (in particular see (p3|)). The driving term fo is of 
non-perturbative origin and is assumed to contribute only for k 2 < q^. The remaining terms 
contribute in the k\ > q$ domain. The driving term thus gives the non-perturbative starting 
gluon distribution 

xg{x,ql) = J° -^-f (x,k 2 ). (5) 

This angular-ordered equation for /, which embodies both DGLAP and BFKL evolution, 
is shown schematically in Fig. 1. The first term in the square brackets in ([|) describes real 
gluon emission with angular ordering imposed. The term containing f(x,k 2 ,q 2 ) is related to 
the virtual corrections corresponding to the unfolded Sudakov form factor, whilst the last term 
in (|4]) represents the virtual corrections which, when resummed, give rise to gluon reggeization. 
The latter correspond to the BFKL part of the (unfolded) non- Sudakov form factor. Angular 
ordering along the chain, a portion of which is shown in Fig. 2, requires 

Zn-xQn-i < Qn where q n = q tn /{l-z n ) (6) 

and z n = x n /x n -\. The angular ordering continues up to a maximum angle with the limit 
expressed in terms of the hard scale via 0(/z — qz). When z is away from the z ~ and z ~ 1 

2 The folded form (which actually is the CCFM equation [B El) contains Sudakov and non-Sudakov form 
factors, which arise from the resummation of virtual corrections and screen the singularities as z — > 1 and z — > 
respectively. 
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Figure 1: Schematic representation of the CCFM evolution equation, (|j), for the unintegrated gluon 
distribution f(x, k 2 , n 2 ). The variable q is defined by q = qt/(l — z) where z = x/x', and k' t = 
k t + {l-z)q. 



domains, angular ordering is equivalent to the strong kt ordering of pure DGLAP LO evolution. 
At first sight there appear to be three types of large logarithms in (fj). First the usual DGLAP 
logarithms coming from the 

= \k t +(l-z)q\ < q (7) 

domain. Second there are the BFKL-type log(l/x) contributions originating from the 1/z part 
of the real emission term in (|4]) and the gluon reggeization contribution. These two terms can 
be combined together in the function 



11 f(-,\k, + q\ 2 , q 2 ) - B{kl-tf)f(-,e t ,f 



k t + q\ 2 \z J ' \z 



(8) 



where we have assumed z <C 1 , in k' t of (0) . Finally there is a danger that in ([|) we have a large 
logarithm from the region q 2 <C k\ . However, we see from (|4]) that the function f(x, k 2 , q 2 ) only 
extends into the region k 2 > q 2 as a result of the BFKL small x effects, which are subleading 
at finite x. That is, to leading log accuracy, it can be shown, in the so-called "single- loop" 
approximation, that if the last term in (|]) is neglected and 0(/x — qz) is replaced by 0(/a — q), 
then the function f(x,k 2 ,q 2 ) vanishes for k 2 > q 2 . Thus we limit the integration regions in 
(U) to the strongly-ordered domain k' 2 q 2 for those contributions in which the unintegrated 
gluon is multiplied by the part P{z) of the splitting function P{z) which is non-singular at low 
z. That is 

2 /V 

P(z) = P{z) - — * (9) 

It should be noted that the BFKL part, (§), of (for which this approximation is not justified) 
is free from singularity as q — > 0, since the potential singularity of the real emission term is 
cancelled by the virtual contribution. 
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Figure 2: A portion of the evolution chain. Angular ordering requires z n _ig n _i < q n , where q n = 
qtn/0- ~ z n ) and z n = x n /x n -i. 

So finally we have just the large logarithms coming from either k' t <C k t or from 2 C 1. Our 
aim is to develop an approximate treatment of the CCFM equation which incorporates both 
types of large logarithms. 



3 Simplification of the CCFM equation 



To simplify the angular-ordered equation, (H), we rearrange the equation and retain only terms 
which generate large logarithms. To achieve this it is convenient to add and subtract the term 



^ f 1 dz [ 9 ( 1 - z - q A 6(A* - q) Q(k 2 - k?) zP(z) % f(x, K\ q 2 ) 



(10) 



from the right-hand-side of (0), and to group together contributions containing the singular 
2Nc/z part of the splitting function 



P(z) 



P(z) + 



2N, 



c 



(11) 



In this way we obtain the approximate form 

f(x, k 2 , /i 2 ) = f (x, k 2 t ) + g jf * dzQ(k t - q )Q L - 



+ 



x 



X 



?(z)9( Z -x)M^ ^- )-zP(z)xg[x, 



1-z 



2tt jo 
dk' t 2 



dz 



dq 2 



1-z 



, /2 e(k 2 -k?)f(x,kf,q 2 )^S (q 



Q(n-q) zP{z) 



+ F{x,k 2 ), 
(12) 
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where the BFKL-type log(l/x) contribution F(x,k 2 ) is given by (||). The second term on 
the right-hand-side is the pure DGLAP contribution in the large x limit. It comes from the 
strongly-ordered configuration 



k' 2 = \k t + (l-z)q\ 2 < q 2 (13) 

in the second term on the right-hand-side of (|4]). In this configuration the variable q becomes 
k t /(l — z). We have also made the large z approximation such that 



P(z) Q(n - qz) w P{z) 6(/i - q) » F(z) /i 



(14) 



With these approximations we may rewrite the second term using 



/ 



fc ?/(i-*) 2 dk 



12 



- 12 I k-t 



Q(z-x) P(z)f\-,k' 2 { _ 



- zP(z)f\x,k' t 



rl 



k t 



1 - z 



H(-.-.v)P(~)-g[ X [ /'' 



zP(z) xg x, 



1 - ^ 



(15) 



where the upper limit kf of the dk' 2 integration has, to leading log fi 2 accuracy, been replaced 
by - z)) 2 . 

Finally, the third term on the right-hand-side of (|i~2D corresponds to the difference between 
([10| ) and the virtual Sudakov (DGLAP-type) contribution given by the second term in the 
square brackets in (HI), that is to the contribution 



f2 /' */*fe(i-,-*) eo.-,)^) 

27T 7o 7 7r<^ y g y 



^ 2 f(x,k' 2 ,q 2 ) - f(x,k 2 ,q 2 ^ 



(16) 



The first integral is evaluated using the strongly-ordered configuration (k' t 2 k 2 ,q « (1— 2)), 
while the second integration is restricted to the region g(l — 2) > k t . It should be noted that 
for q(l — z) -C &t the two terms in the integrand of ( [TBI ) cancel, since then fc£ fc^. 

The contribution (|i~6|) represents the virtual corrections which have to be resummed. The 
resummation is performed in the Appendix. We obtain 

d 



f(x,k t ,fj, ^ 



9 In k 2 



T g (kt,fi) xg(x, k t 



+ 



dq 2 



dk 12 dL(x,k' 2 ,q 2 ^ 



where 
L(x,k 2 ,n 2 ) 



^Q(k t -q ) f dzQ ( fM— — ^ — 
2tc Jo \ 1 — z 



k' 2 d In q 2 



P{z)e{z _ x) *A*(jL_ 



(17) 



zP(z) xg x, 



1 - z 



(18) 



and where the Sudakov form factor 

T.^^l-tf^J^Vw). (1.) 

The cut-off z' < 1 — kt/p enters on account of the kinematic structure of the real emission term, 
where the upper limit of the p integration is given by k t /(l — z'). Note that T g (q, p) therefore 
implicitly depends on k t . The integration limits defining the Sudakov form factor should be 
understood as arising from the G function constraints 

&(l-k t /p)&(p-p)e(p-q). (20) 

This implies T g (q,p) = 1 if these constraints are not satisfied. In particular T g — 1 for k t > p 
or q > p. 

A nice feature of the result fll7|) is that the unintegrated gluon f(x, k 2 , p 2 ) is entirely specified 
in terms of the integrated gluon xg. The next step is to introduce a single scale unified equation, 
which embodies both BFKL and DGLAP type effects, to determine xg. 



4 Strategy for determining the unintegrated gluon 

We have emphasized that the angular-ordered equation ( |T2"D is a 'two-scale' evolution equation 
for f(x, k^, p 2 ). That is the scales k 2 and p 2 are intertwined by angular ordering. In the previous 
section we have shown how the two-scale unintegrated gluon f(x,k 2 , p 2 ) can be determined once 
we know the integrated gluon xg. Here we describe the procedure to obtain xg from a unified 
evolution equation for a single scale auxiliary distribution 



h(x,p 2 ) = 




If we integrate both sides of equation (|T^) over k 2 up to p 2 , differentiate with respect to logp 2 , 
then we find that h(x,p 2 ) satisfies the evolution equation 

h(x,p 2 ) = g jf 1 "* 7 " dz{e(z-x)P(z)Zg(Z,/) 

-zP(z)xg(x,p 2 )} + F(x,k 2 = p 2 ), (22) 

where F(x,k 2 ) is defined by (§). Note that the integral in ( ^2]) has no singularity close to 
z = 1. We can now derive a relation expressing the two-scale unintegrated gluon distribution 
f(x,k 2 ,p 2 ) in terms of the one-scale distribution xg (or h). From (|T7|) we get the following 
expression for f(x,k 2 ,p): 

f(x,k 2 ,p 2 ) = T g (k u p)h(x,k 2 ) + T g {k u p) r^'dzzPiz) ggfflZS - Z)2) xg ( x ,k 2 ) 

JO 271 
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a s rl—qo/kt 

^T g (k t ,n)] o dz 



Q {z _ X )p( z )lg g, kfj - P(z)zxg (x, k 2 ) 



l-kt/n as 



dz — — T, 



2n' 9 \l-z 1 



G(z-x)P(z)-g(- ' kl 



— P(z)zxg x, 



z \ z \ 1 — z 

2^ 



1 - Z 



(23) 



where the scale of as is taken to be the scale of the appropriate gluon, except for the second 
term on the right-hand-side. We may safely set the 1 — z cut-off qo/k t to zero, but not the 
cut-offs k t j 1 \x. 

Note that in the leading ln(l/x) approximation we may set T g = 1 and neglect all the 
integral terms in (^3|), since they do not generate ln(l/x) contributions. In this approximation 
the unintegrated gluon is simply 



f(x, [if 



with no dependence on the scale \i. 

When obtaining (p3|) from flT7|), we have chosen to neglect a contribution coming from the 
derivative of the Sudakov form factor, T g (q,fj,), with respect to k t , which arises from the k t 
dependence of the regulator, see (jl9|). For this reason the unintegrated gluon of (^3|) does 
not precisely integrate to xg(x,fi 2 ), although the corrections are subleading in log/z 2 . The 
discrepancy is indeed negligible at low x, but can become of the order of 20% or so for large 
values of x > 0.1. Rather than complicating (^) by including the derivative, we eliminate the 
discrepancy by changing the regulator in the form factor ([19]) from k t /p to q/p, that is we take 
the Sudakov form factor 



r.ta,,) = «p(-yf A 4 a 4f± r**vpM) 

\ Jg 2 p z 2tt Jo J 



(24) 



This approximation is justified since in our case either q = k t or q ~ k t . Within this approxi- 
mation it is evident that the unintegrated gluon (|23|) integrates exactly to xg(x,fi 2 ); the sum 
of the first two terms on the right-hand-side of (|23|) form the total derivative 



d 

dkl 



T g (k t ,[i) xg(x, k 



(25) 



and the integrals of the third and fourth terms cancel each other. 

From ( pT|) we see that the integrated gluon distribution g can be expressed in terms of h, 
namely 

d^l 2 



xg(x,k t ) = xg(x,q ) + 



h(x, yU 



(26) 
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Equations fl22|) and (p3|) , together with (||) and (p6|) , form a system of coupled equations. If 
we substitute / of (|23| ) into the F term on the right-hand-side of (|22|), and take account of 
(|2T>|), then we obtain an integral equation for h. We may solve this equation for the single-scale 
auxiliary distribution h(x,[i 2 ), and then compute the two-scale unintegrated gluon f(x,k 2 ,ii 2 ) 
from (E21. 



It is convenient to simplify the integral equation ( p2|) for h(x,/j, 2 ) using approximations 
which are valid to leading log accuracy. To be precise we simplify the computation of F(x, k 2 ) 
of (||). First instead of allowing the scale q 2 of / to vary, we note that in (H) the dominant 
values of q 2 are such that q 2 ~ k 2 . Moreover we notice that in the strongly ordered domain, 
k' 2 {= \k t + q\ 2 ) <C k\ , that the first term on the right-hand-side of (§) can be simplified using 

il M. f ( x - + n? a 2 \ ~ [* ^ f (- k' 2 k 2 
nq 2 k? f \z Akt + Ql A ) J k? ; U' *' * 

a; /x 



see ([I]). In the remaining contributions to F(x, k 2 ) of (]§]) we can use ( p3]) to approximate / by 
the first term, 

f(x,k 2 ,q 2 ) « /i(x,fc 2 ), (28) 

noting that g 2 ~ /c 2 and T g (kt, kt) = 1. The other terms in (|23|) give only subleading log(l/x) 
contributions to the BFKL kernel. Since this contribution to F(x, k 2 ) goes beyond the strongly- 
ordered part of the kernel, it is also subleading in log/x 2 . After these approximations, eq. fl2"2] ) 
for h(x,fi 2 ) for /i 2 > g 2 , may be written^ 

^ ,1 



/i(x,/i 2 ) = /i (x,/i 2 ) + / M ^ |e(z - x) P(zU g 2> ) - ^(2) /i(x, g 2 

27r 7o Jql q 1 V \Z J 




where k' 2 = \k t + q\ 2 with /c 2 = /z 2 . The driving term, which arises from the substitution of 
(|26|) for xg in (|22"D and (0), is given by 

/^o^/i 2 ) = £dzi [ e(z-x)P(z)- z g(- z ,q 2 ^ - zP(z) xg(x,fi)}. (30) 

Notice that the strongly-ordered contribution of P( ) has combined with the residual 

DGLAP contribution in (||) with the effect that P{z) -> P(z). 

Equation (23) is the single-scale unified BFKL/DGLAP equation for the gluon that was 
proposed in Ref. [H|. There it was shown that it is straightforward to incorporate a major 

3 In order to be consistent with (p2[), the upper limit of the z integration in the second term on the right- 
hand-side of ( p9| ) should be 1 — qo/q rather than 1. The integrals are, of course, regular at z = 1 and so the 
discrepancy is subleading in In \j? . 
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part of the subleading order log(l/x) (or BFKL) effects by imposing a consistency condition 
to ensure that the virtuality of the exchanged gluon is dominated by its transverse momentum 
squared. This is achieved by including the theta function 6(/i 2 — zq 2 ) in the real emission 
contribution in the last term of (PS|). Other important subleading terms arising from using the 
complete DGLAP splitting function and from the running of as are automatically included in 
our framework Jl4[]. This formalism was used to fit to deep inelastic scattering data was made 
and the auxiliary function h(x,[i 2 ) was determined Ullf . It was checked that the corresponding 
integrated gluon xg(x, k 2 ) computed from (|22| ) was compatible with the gluons obtained in the 
MRS, CTEQ global parton analyses Hl5[ . 

In ]14| the contribution of the quark distributions was included in (|29|) for h(x,fi 2 ). To 



incorporate the quarks in the present analysis we must also include the contribution of the 
singlet quark distribution £ in (|23| ) for f(x, k 2 ,fi 2 ). That is we make the replacement 



P(z) - g - P(z) - g ( -,// 2 ) + P„(z) ~ £ (J,// 2 ) (31 ) 



in the real emission part of the third and fourth terms on the right-hand-side of (|2"B|), where 
// is the appropriate scale. Recall that P{z) = P gg (z). In the second term, and in the virtual 
part of the third and fourth terms, we make the replacement 

zP(z)xg -> (zP(z) + 2n f zP qg (z)) xg(x^' 2 ), (32) 

where rif is the number of active flavours, and the scale // 2 = k 2 or k 2 /(l — z) 2 as appropriate. 
Finally we have to modify the Sudakov form factor so that ( p4|) becomes 



\ Jq 2 2n p 2 Jo 



dz' . (33) 



The above procedure allows the determination of the approximate solution f(x, k 2 , /i 2 ) of 
the CCFM equation, which incorporates both a full (or so-called "all- loop" |6|) resummation 
of leading ln(l/x) contributions, as well as the resummation of leading ln/x 2 contribution, and 
the inclusion of dominant subleading ln(l/x) terms. 



5 The pure DGLAP limit 

It is informative to compare the predictions above for the unintegrated gluon f(x, k 2 ,fi 2 ) with 
those obtained in the DGLAP (or so-called "single- loop" 0, |(|) approximation, in which 0(/U — 
qz) in (H) is replaced by Q(/j, — q), and the last term in (|j) is neglected. After making these 
modifications we repeat the procedures of Sections 3 and 4 and obtain the DGLAP form 

Jx 271 
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r h W) os(t ' /(1 - 2>2: 

h 2n 



T g (k t ,n)xg{x, k 2 ) -T g 



xg \x, 



1 - z 



(34) 



with / = if kf > //. Apart from the last term, this is the equation for the unintegrated gluon 
introduced in Ref. @. The last term, which is only non-zero on account of different scales, 
introduces subleading corrections. Its inclusion improves the accuracy of the integration of 
f(x, k 2 ,[i 2 ) to reproduce xg(x,fi 2 ), see ([I]). Note that the DGLAP or 'single-loop' unintegrated 
gluon vanishes for k t > fi, as indeed it must [§, ||. In the results presented below we include 
the quark contributions as described in (pl|)- (|33|) . 

It is informative to see how the full equation, (123]), for /(re, /c 2 , //) reduces to the DGLAP 



limit (|3^). A crucial observation is that in the DGLAP domain (k t < fi) it is possible, within 
the leading ln(l/x) and ln(/x 2 ) accuracy, to replace P{z) by P{z) in (|23"D . Thus, provided k t < yU, 
we obtain the following more symmetric formula 

f(x,kln 2 ) = T 9 (k t ,ti)h(x,k 2 ) + T g (k t ,») r kt/ "dzzP(z) ^MI£ ~ *)*) xgixtf) 

JO Z71 



a s (k 



2tt 



l—qo/kt 



IJ T g (k t ,n) I dzP(z) 



Q(z- xy-g (^-, kfj - zxg (x, kf) 



+ Q(ji-kt) I dz^Tg 

Z7T 



k t 



1-z 
zxg I x, 



,li)P{z) 



®{z-x)-g ' 

z \ z \ 1 — z 



k t 



(35) 



In the DGLAP limit the first and third terms on the right-hand-side of (^) exactly cancel^; 



they simply represent the DGLAP equation for h(x, kq) of (£1]). Thus (^) reduces to (|34|). 



6 Numerical evaluation of the unintegrated gluon 

In Fig. 3 we show the k t distributions of three different unintegrated gluons at each of four 
different values of x at a hard scale /j 2 = 100 GeV 2 . Briefly 

(i) the continuous curves are the gluons f(x, k 2 , ft 2 ) obtained from (|3|), with the quark 
terms included, using the auxiliary function h(x,k 2 ) of Ref. j[4j], which was itself 
obtained from a fit to deep inelastic scattering data using a unified BFKL/DGLAP 
equation. (The curves have been smoothed in the transition region k t ~ yU.) 

4 Strictly speaking the cancellation is only exact when q — > in (|35|). 
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(ii) the dot-dashed curves show h(x, k 2 ) itself which is independent of /i 2 , 



(hi) the dashed curves show f(x,k 2 ,fi 2 ) calculated from the pure DGLAP equation, 
(|35|), using in this case the auxiliary function h(x,k 2 ) obtained in |14| from pure 
DGLAP evolution from exactly the same starting distributions (xg(x, g 2 ) etc) as 
those found in the unified fit. 
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Figure 3: The continuous curves show the kt dependence of the unintegrated gluon distribution 
f{x,klfi 2 ) for fi 2 = 100 GeV 2 . For comparison we also show the input auxiliary function h(x,k 2 ) 



(dot-dashed curves) [14| and the kt dependence coming from pure DGLAP evolution (dashed curves). 



Note that the third set of gluons are shown solely to illustrate the difference between two 
types of evolution. The gluons of the third set have not been constrained by a fit to the data, 
so should not be regarded as realistic. 
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In the pure DGLAP case, (iii), we see that the distributions are confined to the domain 
k t < /i, as anticipated from strong ordering. On the other hand the distributions f(x,k 2 ,/i 2 ) 
obtained in the unified BFKL/DGLAP framework develop a more and more extensive k t > fx 
tail as x decreases. At small k t and low x the magnitude of the unintegrated gluon calculated 
from the unified scheme is about a factor of two less than that of the gluon coming from the 
pure DGLAP approach of case (iii). This is due to the imposition of the consistency constraint 
in case (i) which suppresses the magnitude of the gluon. If this constraint were absent the 
distributions of cases (i) and (iii) would not be that different. We note that the auxiliary 
function h of case (ii) remains different from the unintegrated gluon / of case (i) down to very 
small values of x. 

For kt > n we see that / is greater than h, whereas the DGLAP-driven unintegrated gluon 



vanishes, as it must. In this domain inspection of (22) and (23) shows that / comes purely 
from the BFKL contribution, 

f(x,k 2 ,^ 2 ) = F(x,k 2 ). (36) 

On the other hand h is smaller than / due to the negative contribution of the integral term 
in (|22D . The latter is a DGLAP contribution which is ruled out when angular ordering is 
imposed. It is negative because the 2N c /z contribution has been subtracted from the real 
emission contribution, but not from the virtual term. We see that for x < 0.01 that there is 
about a factor of 2 discrepancy between h and the true unintegrated gluon /. 



7 Summary 

Here we have addressed the issue of obtaining a reliable determination of the (scale dependent) 
gluon distribution f(x, k 2 ,/j, 2 ), unintegrated over the gluon transverse momentum kt, where \x 
denotes the hard scale of the probe. In the leading log(l/x) approximation the distribution is 
given simply by the derivative of the unintegrated gluon with respect to its scale A = k t , see (0), 
and satisfies the BFKL equation. We correct this simple relation by going beyond the leading 
log(l/x) approximation to include both subleading contributions and DGLAP effects. The 
final result for f(x, k 2 ,fi 2 ) is given in (P3"|). To obtain this result we use the appropriate gluon 
cascade formalism based on angular ordering, which leads to the CCFM equation embodying 
both BFKL and DGLAP evolution. It is important to note that the CCFM equation gives 
a well-defined framework to calculate the very quantity that we seek: the unintegrated gluon 
distribution f(x,k 2 ,fi 2 ). Using this formalism we devise a procedure to determine f(x,k 2 ,ii 2 ) 
from the integrated gluon distribution xg(x,q 2 ), its derivative h(x,q 2 ) and the Sudakov form 
factor T g (q,p,), cf. (p3|). An important ingredient is the solution of a (single scale) evolution 
equation for h(x, q 2 ) which embodies both BFKL and DGLAP effects. From the low x viewpoint 
it includes subleading effects from (i) the consistency constraint which limits the available phase 
space to the region in which the virtuality of the exchanged gluon is dominated by its transverse 
momentum squared, (ii) DGLAP effects generated by that part of the splitting function P gg (z) 
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which is not singular in the limit 2: — *> 0, (iii) the inclusion of the quark contribution and, (iv) 
allowing the coupling as to run and depend on the local scale(s) characteristic of the vertices 
of the cascade. 

We presented sample results to show that the structure of the kt distribution of the gluon 
f(x,k 2 ,fi 2 ) can be significantly different from that of h(x,k 2 ), down to very small values of 
x. There are important consequences for the description of hadron-initiated hard processes in 
which the k t of the gluon is probed locally. 



Appendix 

Eq. fll7D for the unintegrated gluon f(x, k 2 , /i 2 ) is obtained from (|12|) by resumming the virtual 
corrections given in the third term on the right-hand-side of fll2|) . Here we show how the 
resummation is performed. First we note that this virtual correction term can be written as a 
derivative, that is 



d ^Q{kt-k?)f{xXA 2 )\s (q 



f(x,k 2 ,q 2 )Q[q 



-k 2 t 



d 

dk 2 



1 - z 



A{k t , q) R{x, k 2 t , q 2 ' 



k t 



where 



R(x,k 2 tJ q 2 ) 



k\ dk 12 , 2 

-j-pj- f(x, k t ,q ), 



A(k t ,q) = ^ /' dzzP{z)e(l-z-k t /q). 

Z7T JO 



(37) 



(38) 



(39) 



The function R(x, k 2 , q ) has a simple physical meaning. It is the gluon distribution for fixed 
impact parameter b ~ l/k t at scale q. Note that if b = 1/q then the distribution R reduces 
to the integrated gluon xg (x, q 2 ). Using (|37|) we see that ( |12D can be expressed as an integro- 
differential equation for R 



f(x,k t ,fi ) 



2 dR(x,k 2 ,fi 2 ) 
H dk 2 



L(x,^,/i 2 ) + F(x,kf) - k 



, d [» 2 dq 2 



A(k t ,q) R(x, ki,q 2 ), 



(40) 



dk 2 Jh» q z 

where L and F are defined by ( |TH| ) and (|^) respectively. 

In order to solve (|40D for R we integrate both sides over dk 2 /k 2 up to k 2 and obtain the 
following integral equation 



R(x, k t ,n) 



xg (x) + 



*? dk' 2 



il K 2 
» 2 df 

fc? q 2 
14 



L(x,k' 2 ^ 2 ) + F(x,k 



12- 



A(k t ,q) R(x,k 2 ,q 2 ). 



(41) 



From (21) and (22) we see, if \x = k t , that the first two terms on the right-hand-side of (|¥ID are 
just xg(x, kf). The solution of ( |4"I]) may be therefore written 



R(x,k;,/r) = T g (k t , /j) xg(x,k t ) 



+ £¥ T ^« 2 i{£! f 4 (42) 



where the Sudakov form factor 



T 9 (g,/i) = exp (- 



dp 2 

p2 



(43) 



is in agreement with (|1^). Eq. (|23|) then follows from (E3) after differentiation by d/d\nk 
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